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1. Introduction 

In recent years there has been substantial activity studying the relationships between 
integrable systems, which can be exactly solved via Bethe ansatz methods, and the 
spectra of differential equations [1-4]. See the review article [5] for further references. 
This has lead to a deep understanding of the spectral properties of certain Schrodinger 
operators which are not Hermitian, but possess the more general property of PT- 
symmetry [6, 7]. An initially surprising result was the establishment of reality for 
the energy levels of particular non- Hermitian operators with PT-symmetry [8,9]. An 
alternative approach to prove reality of a PT-symmetric operator is to construct an 
equivalent Hermitian Hamiltonian [10-13]. Such a construction has proven difficult and 
only a few explicit (non-perturbative) results are known [14-20]. For a detailed review 
of the field we refer to [5,21,22]. 

Our goal here is to further expose spectral equivalences between Schrodinger 
operators which are quasi exactly solvable (QES) [23,24], i.e. operators for which 
part of the spectrum can be determined algebraically. We will consider two cases, 
viz. one associated with sextic potentials with an angular momentum like term, and 
another where the potentials are expressed in terms of hyperbolic functions. For both 
cases the starting point is to begin with a Hamiltonian which admits two exact Bethe 
ansatz solutions. Here, the Bethe ansatz solutions are of the Gaudin (additive) form 
with finitely many roots, so our approach is different from the cases [1-4, 8, 9] for 
(multiplicative) Bethe ansatz equations with an infinite number of roots. Each solution 
can be mapped to the QES sector of a one-dimensional Schrodinger operator, and these 
turn out to have different potentials. Since the QES eigenvalues of the Schrodinger 
operators are the same as the Hamiltonian, this establishes a spectral equivalence at the 
level of the QES sectors. 

For the sextic case our results show there is a spectral equivalence of the QES levels 
between certain Hermitian and a PT-symmetric potentials. This provides a starting 
point to determine equivalences in a more general context, which we also discuss. In 
some instances we can establish these equivalences rigorously, using the techniques of 
Bender-Dunne polynomials [25] and Darboux-Crum [26,27] transformations. In other 
cases we give conjectures which we find are supported by numerical calculation of the 
spectra. 

In the last section of the paper we study a second class of Hamiltonians from which 
we can determine spectral equivalences of the QES levels for hyperbolic potentials. Here 
we find that the equivalence exists between two Hermitian potentials or alternatively, by 
use of a unitary transformation, two PT-symmetric potentials. This contrasts the case 
of the sextic potentials where the equivalence is between Hermitian and PT-symmetric 
potentials. For the normalisable cases, we obtain a complete spectral equivalence by 
showing that the two potentials are supersymmetric partners [28,29]. 
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2. Spectral equivalences in sextic potentials 

We begin by considering a class of Haniiltonians given by 

H, = A,s; + B.idls; + ^pd,s;), (1) 

where {dp, dj^} are boson operators and {Sp, Sp, S'} are su{2) operators, while ^p = ±1 
is a discrete variable. Through use of the algebraic Bethe ansatz (ABA) method in the 
quasi- classical limit (see for example [30]), the energies and corresponding Bethe ansatz 
equations (BAE) are found to be 

Alp 

Ep = Ap{Mp + Kp) + ipBp Y, 4'\ (2) 

Up P j^^j Up Up 

The parameter Hp is determined by the reference state |0p) which satisfies 

Sp |0p) = f^p |0p) , Sp \(j)p) =0 dp \(j)p) = 0, 
and the eigenstates have the form 

Mp 

2.1. Equivalences of QES sectors 
Next we consider the Hamiltonian 

H = t{na — rife — Tic) + ^{ci^hc + ah^c^) (5) 

where {a, a^ : a = a, 6, c} are canonical boson operators with the usual number 
operators n^ = a^a. It is straightforward to verify that this Hamiltonian commutes 
with the conserved operators A^ = 2na + ui, + Uc, K = ui, — n^. To make a connection 
with ([T]), we make the following two realisations of 5^(2) operators, 



(6) 
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-- a'^c, 


s^ = 
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st = 


--^{no 
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Hereafter 
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set fi 
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(7) 
1. (8) 
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^) = I0(?p)) 



For the reference states we choose 

1 



I0(gp)) 



c 



.t^'^p 



|0) 



which leads to the values 

1^1 = -?i/2, 



n^ = (g2 + l)/2. (9) 

Given a solution for H with the Mi Bethe roots {f i } associated with Hi and gi, we 
need to determine the relationship to the solution of the M2 Bethe roots {^2 } for H2 
with g2- From the form of the eigenstates (jl]) we can deduce the values of the conserved 
operators 

AT = Ml + gi = 2M2 + gs, 

K = Mi-qi = -g2. 

Solving this gives 

gi = M2 + gs, Ml = M2. (10) 

This imposes the restriction Mi < gi. (Note that for Mi > gi the eigenstate (j4]) 
vanishes.) For convenience set g2 = g and M2 = M. If E is the energy corresponding 
to the Hamiltonian ([5]) then through (JTj) we have 



E = El 



{M-q) 



E2 + -{2M + q + 3). 



(11) 



At this point we remark that the energy expression ([2]) and associated Bethe ansatz 
equations ([3]) have precisely the form of the exact solution for the QES Schrodinger 
operator with sextic potential [24]. Explicitly, 



- Mx) + Vp{x)^ljp{x) = 0, 



Vpix) = AEp + x^ - 
for the functions 

ipp{x) 

Qp{x) 



2Ap^pX^ + [ep(4M + AKp + 2) + Aiy + 



(12) 
(2/«p-l/2)(2Kp-3/2) 



x^ 



Ik -i 

X " 2 exp 



M„ 



X 

Y 



Ar,+ 



qpX 



Qp{x) 



n (-' - i'') ■ 

Setting E = — 4(i? + e/2(M — g)) and using the relations (l8ll9|lQ|lip we can write 
-ipp{x) + Vp{x)ijp{x) = E ipp{x) 

Vi{x) = x^ + 6ex^ + x^[2{M - g + 1) + 9e^] 
(M + g + l/2)(M + g + 3/2) 



+ 



X'' 



(13) 



V2{x) = x^ + 6ea;^ + x'^[-4:M - 2(g + 2) + 9e^] - 6e(M + 1) 



+ 



(g + l/2)(g-l/2) 



x^ 



(14) 
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and see that the p = 1 case shares the same QES spectrum as the p = 2 case. 

Before proceeding further some remarks are required. Consider the general 
potential 

V = x^ + 26x' + {6^ + a)x' + ^-^^^ + C (15) 

x^ 

where 5, a and / are real parameters and C is a real constant. Usually we require the 

wavefunction %l){x) to be square integrable along the positive real axis while behaving 

at the origin as ip\x^f) ~ x^^^ . If / > —1/2 this defines a Hermitian problem and so 

we will use the subscript H to indicate this type of radial Schrodinger problem. We 

can alternatively consider PT-symmetric boundary conditions [6,7], denoted by the 

subscript PT. In this case, we require the wavefunction to be square integrable along 

a contour in the complex plane C which starts and ends at |a;| = cxo within the Stokes 

sectors defined as wedges in the complex plane of open angle 7r/4, centered at angles 

— 7r/4 and 7r/4. This contour should also avoid the origin whenever l{l + 1) 7^ 0. In view 

of this we see that the p = 1 case of (fT2|) corresponds to the PT-symmetric case while 

p = 2 is the Hermitian case. 

We mention also that the change of variable e ^ — e is equivalent, up to a unitary 

transformation a^ -^ —a\ to the mapping H — » —H. Hence we also have the equivalence 

that the QES spectrum of Vp{x] e) is the negative of the QES spectrum of Vp{x; — e) for 

both p = 1,2] this is the anti-isospectral duality of [31]. 

2.2. Equivalence beyond the QES spectrum 

In this subsection we will describe how the spectral equivalence obtained in the last 
section between the QES eigenvalues of a PT-symmetric and a Hermitian Schrodinger 
problem is in fact a complete spectral equivalence. 

If we set a = aj = —(4 J -|- 1 -|- 2Z) and C = —26 J in (TT5|) then the results of the 
last section prove that the Schrodinger operators with potentials 

Vn = x^ + 26x^ + {6^ - (4J + 1 + 2/))x' + iillil _ 25 J, 



x"^ 



i)\^^o ~ X' 



l+l 



\/pT = x^ + 25x^ + {5'' + 2J-l- 2l)x^ + ^^ + J)^^ + J+^) (16) 

x^ 

are both quasi-exactly solvable and the J = M + 1 exactly-known eigenvalues coincide. 

The QES spectral equivalence flT6l) at (5 = was discovered and proven [8] via 

the ODE/IM correspondence [1,2]. In fact, the 5th spectral equivalence of [8] makes a 

stronger statement: it says that the full spectrum of the PT-symmetric and Hermitian 

problems flT6|) are isospectral, not just the QES levels. Moreover, the equivalence also 

holds away from the special QES points aj. Numerically, we find this equivalence also 

extends to the case when 5 7^ 0. We conjecture that the most general equivalence is 

that 

Vn = x' + 25x' + {5^ + a)x^ + ^-^^^, ^U^o ~ a;'+^ (17) 

Jb 
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n 


En Hermit ian 


^nPT 





7.17030615 


7.17030616 


1 


19.5220637 


19.5220637 


2 


35.2744653 


35.2744654 


3 


53.7929337 


53.7929339 


4 


74.7062464 


74.7062466 



Table 1. The spectrum of ^ and dH]) with 6 = 0.2, a = 0.31, 1 = 0.54. 



is isospectral to 

^PT = x^ + 25x^ + (5^ - i(a + 6/ + 3))x2 + 



(2/ + 3-«)(2/-l 
16a;^ 



-^ + -(a + l + 2/).(l^ 



The approach used in the previous section does not allow us to prove this statement 
away from the QES sector for the special points aj. It would be interesting to generalise 
the approach of [8] in order to obtain a rigorous proof of this statement. Numerical 
confirmation of the spectral equivalence (IT71) and (IT8l) for values of (5, a, /) away from 
the QES points is shown in table [TJ 



2.3. Further spectral equivalences 

Here we briefly comment that the further spectral equivalences obtained in [8] for ( TT5l) 
with 5 = can also be generalised to the problem with 5 7^ 0. 

The second equivalence of [8] relating the spectrum of two Hermitian sextic 
problems is easily generalised using the results of the last section. The PT-symmetric 
problem flTSl) is invariant under (5, a, /) —* (5, (6/ + 3 — a;)/2, [a + 21 — l)/4) whereas the 
same transformation on the Hermitian case (IT71) leads to a different Hermitian problem. 
Rewriting, we obtain a full spectral equivalence between two radial problems: 



^H^ 



x^ + 25x^ + (52 + a)x^ + 



/(/ + 1) 



-!<^ 



21 -a] 



V'U^o ~ x' 



i+i 



K? = x^ + 25x^ + (5^ + -(3 - a + 60)x^ + ^^^' + "^^'^ + ' + "^ 



'H 



1 

2' 



-^(^ 



16x^ 

(a+2«-l)/4+l 



21 -a), -iplx^o^x 

This equivalence can be proven in terms of Bender-Dunne polynomials [25]. Set 



ipi^x, E, 6, a, I) 



-x^/4:-Sx^/2 l+l 

Jb 



n=0 



-1\" Pn{E,5,a,l) 
a) n\T{n + I + 3/2) 



X 



2n 



(19) 



(20) 



To satisfy a radial ODE with general potential (ITSll the polynomials P„ must satisfy the 
three-term recursion relation 

Pn{E) = {E + C- 6(21 + An- l))Pn-i{E) 

+ 16(n - l)(n + {a + 2l- 3)/4)(n + I - 1/2)P„_2 
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with Pq{E) = 1,Pq{E) = E. The wavefunction ip is an everywhere-convergent series 
for all / 7^ —n — 3/2, n G Z for arbitrary a, with the required behaviour at the origin 
built-in: that is, iplx) ~ x'^^ at x = 0. For generic values of a, the wavefunction i/; 
will be an infinite series. However, whenever a = —(4 J + 21 + 1) for positive integer J, 
the series truncates and the zeros of Pj{E) are the QES eigenvalues [25]. The relevant 
point here is that the recursion relation is invariant under 

a -^ (6/+3-a)/2 , / -^ (a+2/-l)/4 , C -^ C-6{l+2l-a)/2, (21) 

giving rise, with C = —6{1 + 21 — «)/4, to the spectral equivalence ( TT9l) . 

One more spectral equivalence may be obtained, either via the Bender-Dunne 
polynomials as in [8] or using Darboux-Crum transformations [26,27], generalising the 
third spectral equivalence of [8]. Setting aj = —(4 J + 2/ + 1) we take the Hermitian 
QES problem and using a Darboux-Crum transformation remove all of the QES levels 
leaving the rest of the spectrum in place. The resulting potential is once again exactly 
the same as the original form modulo a change in the parameters. That is, the potential 

Vu(x) = x^ + 26x^ + (6^ - UJ + 21 + l))x^ + ^^^ \^\ 

x^ 

V'U^o-a^'+i (22) 



is isospectral to 



Vuix) =x^ + 25x^ + (5^ + (2J -21- l))x^ + (^ + J)(^ + ^ + '^) 

x^ 

+ 25 J, ^U^o ~ x'^-^ (23) 



except for the first J QES levels of (l22l) . The result is unexpected because the 
intermediate potentials, found by removing one energy level at a time, are in general 
singular potentials. Note that this result can also be obtained from the cubic case of 
the type A A/'-fold supersymmetry of [32] . 

Finally, combining the above equivalences, we find that the QES problem 

Vu = x^ + 2Sx'' + iS^-iAJ + 2l + l))x^ + ^-^!-^, V^U^o~a;'+^ (24) 

x^ 

is isospectral to 

VPT = x^ + 26x* + {6^ - (4J + 21 + l))x2 + i^llil (25) 

X'^ 

except for the QES eigenvalues. Numerical confirmation is shown in table [2l 

3. Spectral equivalences in hyperbolic potentials 

Now we move on to examine another case in which a spectral equivalence can be 
estabhshed for potentials expressed in terms of hyperbolic functions. Here we consider 
the following Hamiltonian 

H = e{ni + n2 — n^ — n^) + (7(^1^3 + n2n4 + 01020304 + O1O2O3O4), (26) 
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n 


En Hermit ian 


^nPT 





-11.0798088 


30.1033297 


1 


1.45911359 


48.4085577 


2 


14.4686962 


69.0856538 


3 


30.1033293 


91.8988711 


4 


48.4085576 


116.668373 



Table 2. The spectrum of ^ and ^ with 5 = 0.2, J ^3,1^ 0.54 



where the {oj, a] \j = 1,2,3,4} are canonical boson operators and Uj = a^aj. We note 
that the unitary transformation 



J 



J 



a[ < > Og, a2 < * a^ 



(27) 



is equivalent to the change of variable 

e -^ -e. (28) 

We make the following assignment of the su{2) operators with a central extension: 

S+ = aJ4, S{ = -aitts, 5*1 = -{ui + ns + 1), Ki = -{ui - ns). 



S' 



-a^a^, 02 



0304, 82 = -{us + n^ + l] 



, -2 = 2'^^3-^4j, 



O3 — Ct]^Ct4, 



:+ _.t 



53 = ajai, 53 = -(ni 



"-4 J 



K2 
K3 



-{ni + n^}, 



^4 — 030-2) '-'4 



O2O35 



S'. 



(^3-^2), 



^4 = -(^2 + ^3), 



[K„ S] 



J, -J 



which for j = 1, 2, 3, 4 satisfy the commutation relations: 

I- j ' j -I ~ j ' I- j ' i -I ~ i' I- J' : 
We will also need the corresponding Casimir invariants 

Gj = Oj Oj -\- Oj^Oj I ). 

In terms of these operators we can express the Hamiltonian 

H = H^- g{C, + C2 - {SI + 5|)2 - 27^17^2 + 2{Sl + 5| 



0. 



as 



\-^ 



Hp - g{C, + C4 - (^3^ + SlY - 2K,K, + S^ + ^, 



where 



(29) 
(30) 



H^ = 2e{St - 5|) + g{StS{ + S^S, + 5+^2- + S^S{), 
Hfj = 2e[S^ — 64) + g{S^ 63 + D4 D4 + 63 64 + D4 D3 ). 

We recognise from [30] (see eq. (69)) that Hoi.Hp are Bethe ansatz solvable. Since the 
additional terms appearing in ( l29l) commute with if^, we can extend it to a solution 
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for H. Likewise, since the terms in fl30p commute with Hp, we can also obtain a second 
solution from this expression. 

To obtain the Bethe ansatz solutions we need to identify a reference state |0) which 
satisfies 

SJ 10) = 0, S^ 10) = -s, 10) 

for some scalars Sj] i.e. |0) is a lowest weight state for all realisations of the su{2) 
algebra. The choices 

|0(p.,g.)) = -^^{alr^alr |0) , (31) 

a = a, (3 satisfy this condition with 

Pa + l ga + 1 qp Pf3 

si = —, s, = —, s, = -, s, = -. 

Now we can write down the exact solution from [30]. For the representation (|29ll 
we have that the energies of Ha are given by 

E^ = {p^-q^)e-2Y,v[^^ (32) 

i=i 

where the {va } are solutions of the BAE 



Pa + 1 go + 1 



M 



For a given Bethe root, the eigenstate is given by 

M, 



J2^) W' 3 = l,...,M^. (33) 



t-i \Va + e Va - e/ 



Using this explicit form for the eigenstates we then deduce that 

(5r + 5^) \{v^j^}) = (a4 + 1 + ^^) |{t;i^-)}> 



Paqa I r fi")-. 

Hence the energy of the Hamiltonian ( !29l) in terms of the Bethe roots {va'} is 



i^l^2|{-i^)}> =^|{^i^-)}> 



,0> 



E = ^(p„ + M„)(g„ + M„) + (p„ - g„)e - 2 ^ v^J\ (35) 

For the representation fl30|) the energies of the Hamiltonian if^ are as follows, 

Ep = {pp-qp)e-2Y,vf (36) 
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10 



where the {vn } are solutions of the BAE 



( J 
2 



M« 



+ 



9/3 



,(i) 



+ 



PP 



9 v'^' + e 



,/j) 

-'(3 






(i) 



,W 



J = l, 



M« 



For a given solution, the eigenstate is given by 



{4"}) = n 



S^ 



fc = l 



vy + e 



+ 



b/i 



(k) 
vy -e 



\^iPp,Qf3))- 



(37) 



(38) 



Note that if 

M/3 > min(p/3,g/3), 

then (!38|) vanishes. This is in contrast to (134|) . for which there is no analogous constraint 
on Ma- Taking the above form for the eigenstates, we find that 






qyqp + 2) 



4 



PisiPp + 2) 



K.K. 



{vf}) 



4 
4 



Pl3 + qf3 



{v^^})- 



{V?} 



So that the corresponding energy for the Hamiltonian fl5Ul) in terms of the Bethe roots 
{4^)} is 



Mn 



E = g{Mp - p^)iMfs - qp) - gMp + {pp - q^)e - 2 ^ 



,(i) 



(39) 



i=i 



In order to compare the two Bethe ansatz solutions we need to determine the 
relationship between the parameters {pa, qa, Ma} and {pf^, q^, Mg}. From the form of 
the eigenstates ( l34f38l) it is deduced that 

ni + ni = qa + Ma = qp, 

n2+n3=Pa + Ma= Pl3, 

ni-n2= -pa = Mp-pp 
giving the solution 

qp = qa + Ma, 

PP =Pa + Ma, 
Mp = Ma. 

Hence for every solution of fl33|) with {pa, qa. Ma} giving energy E via fl35l) there is a 
solution of fl37|) with {p^ = Pa + Ma, qp = qa + Ma, Mp = Ma} giving the same energy 
E via ( l39l) . Note that in this correspondence the constraint ( l39l) is never violated. Our 
next goal is to use this result to establish a spectral equivalence in the QES sector of a 
Schrodinger equation. 
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3.1. Equivalences of QES sectors 

We start with the general form of Bethe ansatz equations 

A \~^ ^ //fn\ 

+ v(J) + 7/2 ^ t;0-) - 7/2 " ^ t;0-) - t;(fc) ^ ^ 

for A, 5, C G M and set 

V'(x) = (cosh(s) - l)-(^/2+l/4)^^^gj^^^^ ^ l)-(C/2+l/4) 

M 

(J) 



(4 \ 

— — cosh(x) I TT ( — cosh( 
/ ,=1 



x) + v^'> . (41) 



where 



It can be shown [34] that il){x) satisfies the Schrodinger equation 

- ?T + ^(^)^ = ^V' (42) 

V^(x;A5,C,7) 

= M(M-fi-C+— cosh(s) -l\+-{B + C+lf 

16 4 4 

(25 + 1)(2S + 3) _ (2C + 1)(2C + 3) 

8(cosh(x) - 1) 8(cosh(x) + 1) 

M 

£= -A^v^^\ 

We note that this potential has the symmetry 

V{x- A, B, C, 7) = V{x + iTi; A, C, B, -7). (43) 

Assuming A'y is negative, we see that there is a spectral equivalence between the 
Hermitian problem with potential V{x] A,B,C,'y) and the PT-symmetric problem 
V{x] A, C, B, —7) defined on the contour Q{x) = n. Specifically, there is no PT- 
symmetry breaking in the latter. 

Returning to the Hamiltonian (126|) . we fix (7 = 1. Because the Bethe ansatz 
equations (I33|37l) are of the same form as (HOl) . we can map the spectrum of ( l26l) to that 
of the QES sector of fH2l) by adding the appropriate terms to the potential. Since there 
are two Bethe ansatz solutions for fl26|l we obtain the potentials 

Va{x;pc,, qa) = V{x; 2, -(p« + 1), -(g„ + 1), 2e) + (p« + M„)(g„ + M«) 

+ {Pa - qa)e, 

Vpix; pp, qp) = V{x- 2, qp, pp, 2e) + (Mp - pp){Mp - qp) - Mp 
+ (P/3 - g/3)e, 
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where for Vq,, £ = —2 N^ vj^^ while for Vg, £ = —2 N^ vj^ . We have aheady seen that 

i=i i=i 

the Bethe ansatz solutions for ( 126|) are equivalent when M^ = Mp, pp = Pa + Ma, qp = 

qa + Ma- It follows that the potentials V^(a;;pa, <?«) and Va(x; (p^ + M^), (g^ + Ma)) 

have the same QES spectrum. 

However, the QES wavefunctions of the Schrodinger equation fj42|) are not always 
normalisable on the full real line. The potential ^(x; A, 5, C, 7) has a singularity at 
the origin whenever {2B + 1){2B + 3) 7^ 0. When B = —3/2 or —1/2, the potential is 
nonsingular and the QES wavefunctions (HTj) (assuming A7 is negative) can be extended 
to normalisable odd/even wavefunctions respectively on the full real line [33]. 

It is interesting to consider the two nonsingular cases. The above result with 
{pa, qa) = (—1/2, — 3/2 — M) establishes a QES spectral equivalence between a potential 
with B = —3/2 and a potential with B = —1/2. By adding constant shifts to 
Va(x;pa, qa) and Va(x; (pa + Ma), (qa + Ma)), we can prove that the spectral equivalence 
extends to the full spectrum, except for the presence of a single E = energy level in 
the former. In fact, the potentials are supersymmetric partners [28,29]. Set 

^ , , d (M + l)sinhx . , / \ 



d"^ 

+ V{x; 2, -3/2, M - 1/2, 7) - j{M + 1) 



%l){x) = Eip{x) 



then 

with corresponding E = eigenfunction 

ip(^x) = (coshx + l)(^^+i)/2 ^ g^p ^Zcoshxj . 
We immediately deduce that the supersymmetric partner 
Q_Q+^(x) = E^(x) 

has potential V{x; 2, -1/2, M + 1/2, 7) + M - 7(M + 1). We have therefore established 
complete isospectrality between these Schrodinger problems, up to the E = energy 
level. Moreover, this spectral equivalence holds for all real values of M. 

Finally, we remark that the unitary transformation (l27|l combined with the change 
of variable (128|) leaves the Hamiltonian invariant. Observe that (l27j) applied to (I3T]) has 
the effect that 

p < — > q. (44) 

This is reflected in the symmetry of the Bethe ansatz equations (I33II37P which are 
invariant under the combination of (I28|44j) . Thus the unitary transformation (1271) 
effectively interchanges the Hermitian and PT-symmetric versions of the potentials 
Va{x;pa,qa), V/six] pp, qjs) , which explains the existence of the symmetry fH3|) . 

This scenario is somewhat different to the case of the QES sextic potential discussed 
previously. There, the unitary transformation maps the QES spectrum into the negative 
of the QES spectrum. The equivalence between the Hermitian and PT-symmetric 
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problems is due to the equivalence of the Bethe ansatz solutions. For the above case, 
the unitary transformation maps between the Hermitian and PT-symmetric cases, while 
the equivalence of the Bethe ansatz solutions gives a spectral equivalence between two 
Hermitian QES potentials. 
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